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Introduction. Anomalous baryon and lepton number violating processes in the electroweak 
theory are dominated in the semiclassical weak coupling limit, a w — > 0, by field configu- 
rations which solve the classical Euclidean Euler-Lagrange equations. At zero temperature 
and energy the classical solutions that contribute to anomalous winding number transitions 
are the familiar BPST instantons/anti-instantons [[!]. The rate of these vacuum tunneling 
transitions is exponentially suppressed ~ exp(— 25/), where 25/ = 4n/a w is the Euclidean 
action of a widely separated instanton (I) /anti-instanton (I) pair. At temperatures much 
higher than Myy, the transitions are classical thermal activation transitions over the po- 
tential barrier between vacua, with a Boltzmann rate ~ exp(— i^/fe^T) controlled by the 
energy E s ~ AM w /otw of a certain unstable classical stationary field configuration called 
the sphaleron f2|. 

At intermediate temperatures, or at finite energy (not necessarily arranged in thermal 
equilibrium), very little is known about the rate of anomalous transitions between states of 
different winding number. The first step in studying these transitions is to find the classical 
solutions which dominate the semiclassical rate, and calculate their action. This involves 
solving the classical nonlinear field equations in Euclidean time. At low energies one can 
construct solutions of the classical Euclidean Yang-Mills-Higgs equations consisting of peri- 
odic chains of II pairs arrayed along the imaginary time axis. The action of these periodic 
instanton solutions can be expressed as a power series in (E/E s ) 2/3 for small E/E s ^. This 
perturbative treatment of low energy periodic instanton solutions can be recast as an ex- 
pansion in powers of (My/f3) 2 where f3 is the period of the solution. For larger energies or 
periods the solutions can be found numerically. In this Letter we describe the qualitative 
properties of and present numerical results for these periodic instanton solutions of the SU(2) 
Yang-Mills-Higgs equations, i.e. the bosonic sector of the standard electroweak theory with 
9w = 0. An unexpectedly rich structure of bifurcating periodic instanton solutions, both real 
and complex has been found, whose physical consequences for B and L violating transitions 
in the electroweak theory remains to be more fully investigated. 

Periodic Instantons at Low Energy. As a simple example of a periodic potential with periodic 
instanton solutions consider the pendulum potential, 

V(q) = to 2 (l -cosg) . (1) 

The zero energy instanton which interpolates between the vacuum states at q = and q = 2tt 
is the kink configuration, ?/(r), given by 

cos (^Y^j = -tanh(cur) (2) 

which solves the classical Euclidean equations qj = V'(qi) and has action 5/ = 8u. The 
anti-instanton solution is qj(r) = <?/(— t) with the same action. Consider now the widely 
separated I — I pair configuration, 

qii(r) = qi(r) + qr(r-f)-27r (3) 

with f 3> \ju). The action of this configuration can be computed to first order in the 
interaction between the pair, with the result 
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S[q ir ] = 25/ - 32 uo e' wf + O (V 2 ^) . (4) 

The negative sign reflects the attractive interaction between the I and /. 

We now consider a periodic arrangement of I and I at equal intervals along the imaginary 
time axis, with period (3. This means that the separation between nearest neighbor / and 
I is f = (3/2. The attractive force between neighbors can now exactly balance and yield 
an extremum of the action, the periodic instanton solution. Since there are two nearest 
neighbor II interactions per period we expect the action per period of this solution to be 

S((3) = 16 u - 64 ue~^' 2 + O (e~^) , (5) 

in the limit of large uj(3. Since 

E ^ ) = ^l = 3 2^e"^ 2 + ... (6) 

(XL) 

large (3 corresponds to low energy. In the one dimensional pendulum example the exact 
periodic instanton solution with this action and energy are easily found explicitly in terms 
of elliptic functions by simple quadrature. Because of the attractive interaction between the 
I and I along the chain it is clear that there is a single negative mode of the second order 
fluctuation operator, —d 2 + V"(q(r)), around this periodic instanton solution, a fact that is 
also reflected by the second derivative of the action, 

d 2 S((3) dE{(3) 2 ua/2 

—^r = —nr = - 16uje / +--- <0 - 7 

dp 2 d(3 

The monotonic decrease of period (3 with increasing energy persists up to E = E s = 2a; 2 , 
where the curve of S((3) vs. (3 of the periodic instanton becomes tangent to the constant 
sphaleron solution, corresponding in this simple model with the unstable static configuration 
q g = ii. This occurs aX (3 = (3- = 2ti/uj equal to the period of oscillation in the inverted 
potential at q = q g . At this (3 the action of the periodic instanton is E s (3- = Attuj < 16u, 
reflecting the fact that the action is monotonically decreasing as (3 ranges from oo down to 
and as E increases from to E s . 

Beyond the point where the periodic instanton and sphaleron solutions merge, a complex 
solution bifurcates from the sphaleron and continues with real decreasing action. This may 
be understood by the amplitude of the zero mode at (3 — (3- turning from real to pure 
imaginary as (3 is decreased through the critical value. The generic behavior described here 
is what we call type (I) behavior of the periodic instanton solutions, for which the monotonic 
negative sign in the first derivative of the action in (0) and E s (3- < 2Sj are characteristic. 

A different pattern is possible when the instanton has additional zero modes, and there- 
fore additional parameters enter the description. Such is the case in field theory models 
with exact or softly broken conformal invariance. We have studied this case in some detail 
in the 0(3) nonlinear sigma model in two dimensions, softly broken by a mass term [||. 
This model shares many features with the bosonic sector of the SU(2) electroweak theory. 
Although there is no isolated single I ot I solution in the broken theory, due to Derrick's 
theorem (which tells us that zero scale size p —>■ has minimum action), a periodic instanton 
solution does exist in which the scale size p is adjusted to a certain value as a function of 
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period (3. At this value the attractive interaction between / and I exactly balances the 
tendency of each individual / or / to collapse to zero size. It is again possible to understand 
this at low energies by first finding the two-body interaction between well isolated I and /, 
and then arranging them periodically along the imaginary time axis, calculating the change 
in the action per period from that of a single II pair due to the sum of the first order in- 
teractions between them. In this calculation the scale size p can be treated as a variational 
parameter with the value on the solution p(f3) determined by extremizing S(/3,p) with re- 
spect to p. Substituting into 5* then gives S(/3) on a low energy periodic instanton solution. 
The resultant behavior depends upon the curvature of S(/3), which is of opposite sign relative 
to the type (I) models, and we will consequently refer to this case as type (II), i.e. 

The periodic instanton solutions in this case have two negative modes rather than just 
one, with the second negative mode corresponding to variation of the scale size p away 
from its extremal value p(f3). Such periodic instanton solutions do not contribute to thermal 
winding number processes at finite temperature. However, they can contribute to anomalous 
finite-energy non-thermal transitions. 

Because of the existence of the conformal mode in SU(2) BPST instantons, we would 
expect the Yang-Mills-Higgs theory to behave qualitatively similar to the 0(3) sigma model, 
and to also be of type (II). Indeed for low energies we can show that this is exactly what 
happens. The action for an isolated pure SU(2) II pair with scale size p separated by 
distance f is 

S« = - ^ + O (£) , (9, 

with Si = 8ir 2 /g 2 = 2n/a the single instanton action and the second term the well-known 
dipole-dipole attractive interaction between the / and I aligned in an internal 577(2) direc- 
tion. When the periodic chain of I and I separated by (n + is constructed, this leads to 
the total interaction, 

_ QGttV » 1 _ 47r/27rV\ 

When the SU{2) doublet Higgs field is added to the action it can be solved for at leading 
order in the I or I background and gives a contribution, 



(ii) 



which was first calculated by 't Hooft |J. 

This positive contribution expresses the fact that Derrick's theorem drives the single 
isolated I or I scale size to zero p; however, at finite (3 this is opposed by the dipole-dipole 
interaction S- int , and the variational action 
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has an nontrivial extremum at p = p(f3) when dS/dp = 0, or 
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denoting by x the expansion parameter MwP/(2Tt). This stationarity condition implies that 
the various next-to-leading contributions to the action from both the gauge and Higgs fields 
are all of order x 6 . Hence the periodic instanton action, 
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1 + \x* + O (x 6 ) 



for small Mwf3, and 
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dS{(3) 4M 
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X 
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1 + O (x 2 ) 



(14) 



(15) 



is an increasing function of period f3 for small Myy{3. The second derivative of S(/3) is also 
clearly positive. After a rather elaborate calculation ||, one can evaluate the coefficients of 
the next to leading terms in the perturbative expansion. One thus finds 
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AM, 



l + -x 4 + -x 6 + 0(x 
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W 3 
X 
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2 3 
1 + Ax 2 + O (x 4 ) 



(16) 
(17) 



Although the SU(2) Higgs theory starts out at low energy and small (3 behaving like 
the type (II) case, there is an additional independent parameter in the AD gauge theory, 
namely the quartic Higgs self-coupling A, or equivalently the Higgs mass, Mh- Thus, we 
cannot preclude a more complicated behavior at larger energies, depending on the value of 
Mh/Mw- This is indeed what we have found in our numerical study. 

Periodic Instantons at Finite Energy in the SU (2) -Higgs Theory. We consider the SU(2) 
gauge theory with a doublet Higgs field in AD Euclidean space with the action, 



S 



j 2 J d 4 x j^Tr (F^FfJ) + (D^(D^) + A (V$ 



2 2 \ 2 

9 v \ 



d^A u — dyA^ — i\Ap, A u ] 



(0, 



and Afj, = A a ^a a /2. The corresponding classical Euler-Lagrange equations are 

DpF^ + i(D v &) x$-i$ f x (D„$) = 
2A 



$t$ 



g 2 v 2S 



9 Z 



$ = 0, 



(19) 
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where the covariant derivative acting on the A-field in the adjoint representation is 
D^A U = d^Ay —i [A^,A V ], and the x denotes the outer product of the two spinors. We 
use the standard conventions for the Higgs vacuum expectation value v and the self-coupling 
A in which the W- and and Higgs-masses are M w = ^gv and M H = \/2Xv respectively. 

The spherical Ansatz is given by expressing the gauge field A M and the Higgs field $ in 
terms of six real functions do , a% , a , 7 , u and w of r and r: 



Ai(x,r) 



A:fx,r) 



-a [r,T)a x 

1, , N _ „ A a(r, t) 
- [ai(r, r) a ■ xxH — (cr* 



cr • x^i) H e iifc Xj a k \ 



$(x, r) = V2M W [u(r, r) + iw(r, r) a • x] ( 



(20) 



where £ is an arbitrary unit two-component spinor, x is the unit three-vector in the radial 
spatial direction, and a is a three-vector of Pauli matrices. 
Upon substituting ([ 



S 



4:71 



dr 



dr 







into the action fll8f) one finds @ 



+\{XX + 1)# + Re(«x0 2 ) + 4 r 
where the indices now run over and 1 and 



2M^ 



2^2 ( XX ~ l ) 2 



(21) 



X = a + i[y-lj , x = a -1(7-1), 
> = V2M w (u + iw) , <p = V2 M w (u - iw) , 



I 

Dp(j) = (dp - - a„) 



= (<9 M + - a M ) 



The equations of motion for the reduced theory are 



(22) 



(Ax)x - xDuX 

ixx - 1) + \h> 1 X 



+ - r z 
2 

z 

~2 



1 2A 
-^(^) + ^fe + l) + -^r 2 ( 



(A 
2M[ 



4> = i x4> ■ 



(23) 



Note that the overbar on </>, x an d -A denotes changing z — > — i in the definitions (22) above, 
which is the same as complex conjugation only if the six fields 0^,0,7,^ and w are real. 
These equations can be obtained by either imposing the the spherical Ansatz (20) on the 
four dimensional equations fll9"|), or by varying the action (|2~T|) directly. 
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The spherical Ansatz ( ^0|) has a residual U(l) gauge invariance under the U(l) gauge 
transformation, 

X^^X, 

0^e iQ/2 0, (24) 

The complex scalar fields x an d have U(l) charges of 1 and 1/2 respectively, a M is the U(l) 
gauge field, f^ u is the field strength, and is the covariant derivative. The residual U(l) 
gauge invariance must be fixed for numerical solution of the equations. In our numerical 
work we chose the temporal gauge = 0. The remaining time independent gauge freedom 
is fixed by a boundary condition at the quarter period time slice r = (3/4 to be specified 
below. In the a = gauge the v = component of the first of Eqs. (|23D , i.e. the Gauss law 
constraint, must be imposed on the initial r = surface, whereupon it will be satisfied for 
all r. 

The action of the various discrete symmetries, C, P, and T on the two dimensional fields 



follows directly from the spherical Ansatz (p0|). In addition to these symmetries, we may 



consider the two dimensional reflection symmetry R : — > — 0. We can employ these discrete 
symmetries to help us select the appropriate boundary conditions for the periodic instanton 
solution. The CPR even fields are 7 and w, whereas the other four fields are CPR odd. 
Since we are searching for a periodic solution which returns to itself with period (3, the time 
derivatives must reverse sign in the second half period relative to the first. This means 
that we should require the boundary condition that the r derivatives of all remaining five 
functions in gauge vanish at r = and r = (3/2. This corresponds to an instanton at (3/4 
and an anti-instanton at 3(3/4 in the low energy limit. At the time slice r = (3/4 the fields 
are sphaleron-like. In a gauge where the four dimensional fields are regular at the origin, the 
sphaleron is a CPR even configuration, and therefore the CPR odd fields change sign while 
the CPR even fields reach a maximum at r = (3/4. Hence we actually require the solution 
only on the quarter interval [0,/3/4], if we specify the boundary conditions, 

a 1 = a = A f = u = w = 0, r = ; 

(3 , , 

ai = a = u = = j = w, r= 4' ( ) 

These boundary conditions eliminate the time translational zero mode. 



The Euler-Lagrange Eqs.^3] are also invariant under the two additional complex discrete 
transformations: 

d : 0^-0, 0^-0, X^~X, X^-X, a->-a (26) 
C 2 : 0^0*, 0^0*, X - X* , X - X* , a-^a* . (27) 

Under C\. S — > S, under C 2 : S — ► S* (and similarly for E). 

For the four dimensional fields to be regular at the origin, and to approach the vacuum 
at r = 00 we require the boundary conditions, 



a = 7 = w = = a! x — u' , r = 
a\ = a = u = 0, but 7 = 2, w — 1 , r = 00 , (2£ 
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where the last condition is necessary for a nonzero winding number, and agrees with the 
sphaleron boundary condition at r = oo on the r = /3/4 slice. The boundary conditions on 
di at the origin and infinity are gauge choices, which completely eliminate the time indepen- 
dent gauge freedom in temporal gauge a = 0. 

Numerical Results. With a well-defined elliptic boundary value problem standard numerical 
methods may be applied. Here we only outline our computational procedure, the full details 
of which will be presented in a separate publication ||. We discretized the Euler-Lagrange 
equations in space and time by using link variables for the gauge degrees of freedom, thus 
preserving exact gauge invariance under space dependent gauge transformations (consistent 
with our choice of the ao = gauge). The code allows for non- uniform grids in both space 
and time to better approximate the fields in the regions of fastest variation. We also allowed 
for the analytic continuation of the solutions into complex valued functions by using complex 
variables to represent the fields ai,a,j, u, w and by discretizing the equations in a manner 
compatible with analytic continuation. The coordinates (r, r) are maintained real. 

The equations were solved by the Newton-Raphson technique. Denoting a definite trial 
configuration of the fields by fi, where i stands for the discretized two dimensional lattice 
point (r, r) as well as the various field components themselves, we may calculate the gauge 
invariant discretized action functional S[f], the discretized first variation, dS/dfi, and the 
second order fluctuation operator d 2 S/dfidfi' on the trial configuration. Since the first vari- 
ation must vanish on the solution, one can find the first order correction to the configuration, 
Sf, by solving the linear equations, 

Vik.^ + wr - (29) 

Adding 5f to / yields a corrected trial configuration, and the process may now be iterated. 
Clearly, all gauge and translational zero modes must be removed from the second order 
variation by the gauge fixing and boundary conditions in order for the inverse of d 2 S/dfidfi> 
to exist and the procedure to be well defined. 

The algorithm converges quite rapidly, the error decreasing quadratically with the number 
of complete Newton-Raphson iterations. The most time and memory consuming step is the 
inversion of the second order fluctuation operator. With grids consisting of as many as 
128 x 128 points (or more) and 5 complex valued fields per point, a direct solution of the 
above system of linear equations is prohibitive. Because the equations of motion are local, 
the matrix to be inverted is a sparse banded matrix, and it is much more efficient to use 
the method of eliminating alternate time slices instead of a direct inversion. This effectively 
reduces the dimensionality of the linear system one must invert to the size of the space grid 
only, times the 5 components of the fields. Quite modest lattices (64 x 64) are sufficient for 
accuracy of order one percent, except in the small f3, low energy region, where the larger, 
adaptive lattices (128 x 128 or greater) were used. From this solution at a given (3 and A 
other solutions were found by using the previous one as a trial configuration for the new 
values of the parameters, changing the values of period and A in small increments. 

Our results for the dependence of the action on the period are shown in Figs. 1 for two 
different values of the quartic Higgs coupling, A = 0.7g 2 and A = 3.6g 2 , respectively. They 
clearly exhibit a pattern of bifurcations. 
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FIG. 1. The action of the periodic instanton as a function of period along the four distinct 
branches of solutions. The action of the sphaleron is shown by the solid diagonal lines. In (a), 
A = 0-7g 2 < A cr it- The only real solutions lie along branch 1, which joins onto the perturbative 
solutions at small period. Branch 1 merges with the sphaleron at s, and may be analytically 
continued onto branch 2. Branches 3 and 2 have real action and energy and form a cusp at point 
c, beyond which along branch 4 the action becomes complex (only the real part of the action is 
shown). The cusp c always lies above the sphaleron line in this regime of Higgs coupling. In (b), 
A = 3.6g 2 > A cr it and the situation is similar, except that it is now the complex branch 3 that 
merges with the sphaleron. Branches 1 and 2 have real action and energy and form a cusp c that 
lies below the sphaleron line. Again, the action and energy along branch 4 are complex. 

Since the numerical method and our code works equally well for real or complex solutions, 
we were able to follow the action and energy of the latter as well. We observe an interesting 
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pattern of bifurcating solutions depending on the Higgs self-coupling A, as anticipated in ref. 
I- 

There are two regimes, separated by a critical value A cr . As illustrated in Fig. la, for 
A < A cr the perturbative solutions (branch 1) merge with the sphaleron at some value 
[3 = j3\, which is similar in behavior to the 0(3) sigma model. For (3 > f3\, branch 1 can be 
analytically continued up to a second bifurcation point j3 2 to yield complex solutions whose 
energy and action however remain real (branch 2). At 02, branch 2 merges with yet another 
branch of complex solutions (branch 3), and the action and energy of these solutions are also 
real. The numerical results for this second complex branch suggest that its action decreases 
monotonically with decreasing /3, and that at a sufficiently low (but positive) (3 its action 
vanishes, with potentially interesting physical consequences (however, as indicated below, 
the energy along this branch becomes arbitrarily large). Branches 2 and 3 form a cusp c at 
/?2> beyond which the solutions may be analytically continued onto a fourth branch. The 
action and the energy on branch 4 (f3 > /3 2 ) are complex, and the cusp c always lies above 
the sphaleron line (only the real part of the action has been graphed). 

Figure lb illustrates the situation for A > A cr . The pattern is similar, except in this 
case it is the complex solutions along branch 3, and not the real solutions along branch 1, 
that merge with the sphaleron, and whose extension beyond fti becomes branch 2. Branch 1 
(the perturbative one) and branch 2 are real up to /?2, where they form a cusp c beyond 
which the solutions are complex with complex energy and action (branch 4). In this case, 
the cusp always lies below the sphaleron line, except at a critical coupling A cr i t where the 
cusp intersects the sphaleron. We found A cr ~ 1.198g 2 , corresponding to Mh — 3.096Mw, in 
good agreement with 3.091Mjy, obtained in ref. by a careful treatment of perturbations 
away from the static sphaleron solution. 
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FIG. 2. The energy as a function of period near the bifurcation points for two values of A on 
either side of the critical Higgs coupling and for the critical coupling A cr it itself. The solid horizontal 
lines labeled by E s denote the constant sphaleron energy for each value of A, while the dashed, 
dotted, and solid curves are the energies of the corresponding solutions in Figs. 1. 
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The bifurcation structure for the energy of the complex solutions is shown in Fig. |2] for 
two values of A on either side of the critical value and for the critical coupling A cr i t itself. The 
associated sphaleron energies E S (X) are illustrated by the short horizontal lines. The real 
solutions are indicated by the bold solid lines, the complex solutions with real energy by the 
dashed lines, and the complex solutions with complex energy by the dotted lines extending 
down from the cusps c (only the real part of the energy has been graphed). The critical 
coupling A cr it occurs when the energy E (X) of the cusp crosses the sphaleron line, i.e. when 
-^o (Acrit) = E s (Acrit). Finally, it should be noted that the equations are analytic (both in 
A and in (3), and therefore the total number of solutions cannot change at the bifurcation 
points. Indeed, branches 1, 2 and 3 consist of two independent solutions each (with the 
same energy and action), whereas branch 4 consists of 4 solutions (in two conjugate pairs). 
By continuity in A, these (unsuspected) new complex instanton solutions continue to exist 
even at lower A. The physical consequences of these new solutions for rates of anomalous 
processes at finite energy in the electroweak theory are currently under investigation. 
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FIG. 3. The action (a) and energy (b) of the periodic instanton as a function of (3 in the pertur- 
bative regime (small (3). The numerical data points are denoted by the triangles, the leading order 
calculations (14)-(15) by the dashed lines, and the next to leading order calculations (16)-(17) by 
the solid lines. 
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To conclude, let us remark that the accuracy of our numerical calculations is high enough 
that they can provide a verification of the perturbative results for low energy, low period. In 
Fig. we plot our numerical results for action and energy together with the results of the 
perturbative expansions to leading order (|14])-(|15l) and next to leading order ( (L6f )-(|17|). The 
agreement with the next to leading order perturbative results is quite good. 
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